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ABSTRACT 

We investigate the evolution of young star clusters using iV-body simulations. We 
confirm that subvirial and fractal-structured clusters will dynamically mass segregate 
on a short timescale (within 0.5 Myr). We adopt a modified minimum-spanning-tree 
(MST) method to measure the degree of mass segregation, demonstrating that the stars 
escaping from a cluster's potential are important for the temporal dependence of mass 
segregation in the cluster. The form of the initial velocity distribution will also affect 
the degree of mass segregation. If it depends on radius, the outer parts of the cluster 
would expand without undergoing collapse. In velocity space, we find 'inverse mass 
segregation,' which indicates that massive stars have higher velocity dispersions than 
their lower- mass counterparts. 

Subject headings: methods: numerical - open clusters and associations: general - stars: 
kinematics and dynamics 



1. Introduction 



Most stars with masses m > 0.5Mq are thought to form in star clusters (jLada &l Ladall2003l ). 
It is, therefore, very important to understand star formation in the context of star cluster formation. 
Observatio ns of very young and forming star clusters imply dynamically cool and clumpy initial 
conditions ([Williams et alj|2000i ; ICarpenter Hodappl 120081 ) . 



Young clusters (< 6 Myr) also often exhibit 'mass segregation,' where the massive stars are 
more centrally concentrated compared to their lower-mass companions. The origin of this mass seg- 
regation has been suggested as either 'primordi al,' that is, it is a result of the star-formation process 
in which stars (particularly t he massive stars: iMoeckel fc Batd 120101 ) form mass-segregated from 
their parent molecular cloud (jChen et alj 120071 : IWeidner et al.ll2010l ). or dynamical, i.e., resulting 
from fast dynamical evolution, with increased importance in the denser cores. 
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Allison et ajj (|2009al ) investigated the evolution of initially subvirial and highly substructured 
open-cluster-like objects. They found that this substructure will disappear on very sh ort timescales 



(a few Myr), with clusters undergoing mass segregation. Subsequently, these authors ([Allison et al 



20101 ) also found that the more clumpy and cooler the cluster is, the higher the degree of mass 
segregation will be. However, not only the spatial distribution but also the stellar velocity distri- 
bution affects the degree of ma ss segregation: initia l ly sub virial clusters will exhibit gradients in 
their radial velocity structures. iProszkow fc Adama (|2009i ) found the initial subvirial state of the 
star-forming clumps also plays an important role in the sense that it facilitates high interaction 
rates, which leads to rapid mass segregation. 

In this paper, we investigate the influence of different initial velocity distributions on the 
development of mass segregation. In §2, we introduce our initial conditions. In §3, we analyze our 
simulations, and in S4, we discuss the results and draw our conclusions. 



Method and initial conditions 



2.1. Initial cluster model 



Both observations and theory suggest that young star clusters form with subvirial kinematics 



and highly substructured in space and velocity ( I 


parson 


1995: 


Williams et al. 


2000: 


Elmegreen 


2000; 


Testi et al. 


2000; 
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; Carpenter & HodaDD 
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Proszkow & Adamsll2009). We simulated the evolution of cool, clumpv star clusters containing 100C 



single stars and no gas wit h initial virial r adii r v = 0.5 pc. The stellar initial mass function follows a 



three-part power law after iKroupal (120021). with minimum and maxim um masses of 0.08 and 5OM0, 
respectively. We follow the method of lGoodwin &: Whitworth! ([20041) to generate 'fractal' clusters. 
First, we define a cube with sides N^ v . We place the first-generation 'parent' in the center of the 
cube. Next, we split the cube into N^ iv parts, each containing a subnode ('child') in its center. We 
adopt Adiv = 2 throughout, so that there will be eight subcubes and eight children for each parent 
cube. 

A child has a probability of Njj? 3 ^ to become a parent of the next generation, where D is 
the fractal dimension. The lower D is, the fewer children become next-generation parents, and 
the clumpier the cluster will be {D = 3.0 corresponds to a uniform distribution). After removing 
the first-generation parents and the children that do not mature, we add Poissonian noise to the 
positions of the parents to prevent the cluster from developing an obvious and artificial gridded 
structure. Each mature child then divides into A^ iv children in the centers of A^ iv subsubcubes. 
This process is repeated until there are many more children than required. Eventually, cluster 
stars are selected randomly from the remaining children. In this paper, we will not discuss effects 
associated with varying the fract al dimension, so we adopt a fixed fractal dim ension D = 2.0 to 



satisfy observational constraints (jMontuori et al 



1997 



Sanchez & Alfaroll2009l ). 
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Goodwin fc Whitworthl (|2004l ) also suggest that the velocity dispersion must be coherent, which 
means that the children inherit their parents' velocities. A random velocity component is also added 
to each child. We argue that for star formation in molecular clouds, the velocities of the stars depend 
on the potential they are embedded i n, in the sense that the stars in the inner parts should have 
higher velocities. iKiipper et al.l (|2010l ) show that the radial velocity-dispersion profiles they derived 
based on a set of iV-body simulations decline as a function of radius for all bound stars as well as 
for the cluster member stars within the Jacobi radius (where the internal gravitational acceleration 
equals the tidal acceleration from the host galaxy). 



On the other hand. lRochau et al,l (|2010l ) use relative proper motions to show that in the young 
(~ 1 Myr) Galactic starburst cluster NGC 3603, stars with masses between 1.7 and 9M all exhibit 
the same velocity dispersion, which implies that the cluster stars have not yet reached equipartition 
of kinetic energy. Inspired by these constraints, we will therefore consider two different forms of the 
radial velocity-dispersion profile in our iV-body calculations, i.e., a constant and a radially declining 
functional form. Although molecular clouds do not possess a single dense core, there still exists 
a dynamical center surrounded by many molecular-cloud clumps. Thus, we define our radii with 
respect to this dynamical center. The velocity distribution is then written as 



v 



parent 



a = constant or a 



(7 



(1) 

1 + (r/r v p (2) 

where f pa rent is the velocity of the parent node and a the velocity dispersion. Since the radial 
form of the nonconstant initial velocity distribution is not well understood, we adopt a plausible 
functionality that leads to a large velocity dispersion, which in turn facilitates our exploration of 
its effects. We will discuss both kinds of velocity dispersion in §3. 

We only follow the evolution of our simulated clusters for a very short time (< 6 Myr), so that 
no stellar evoluti on is included , given that i ts effect would be negligible on these timescales. We 
assume, following lAllison et al.l (|2009al . |2010| ) and for computational reasons, that our simulations 
commence at the time that all stars have formed. We do not include the effects of the surrounding 
gas to reduce adding to the complexity of the physics involved; we aim at doing so in a subsequent 
paper. 

Finally, for clarity we summarize the other parameters used for our fractal model: 



1. Total number of stars, iV = 1000; all stars are treated as single stars; 

2. Fractal dimension, D = 2.0; 

3. Initial virial ratio, q = 0.30 (here defined as the ratio of kinetic to potential energy); 

4. Virial radius, r v = 0.5 pc. 



- 4 - 



We generate 100 samples for every velocity-distribution model. These samples only vary with 
the initial random seed. This way, we can address t he resulting cluster properties statistically. We 
used kira in starlab (jPortegies Zwart et al.ll200ll ) as our iV-body integrator. 



2.2. Measuring mass segregation 



We use the minimum-spanning-tree (MST) method ([Allison et al.ll2009bl ) to quantify the degree 
of mass segregation. An MST connects all sample points using the shortest path length, without 
forming any closed loops. There can be multiple shortest paths, but the length is unique. We 
measure the degree of mass segregation of the N most massive stars by comparing their MST 
length with the average MST length of N randomly selected stars in the cluster (where both values 
of N are the same). 



Allison et al.l (|2009bl ) define the amount of mass segregation by the ratio of the random to the 



massive-star MST lengths, 

^ (/norm) _j_ ^norm 

/massive /massive 

where A is the level of mass segregation (the 'mass-segregation ratio'), /massive is the MST length of 
the massive stars, and (/ n0 rm) an d cr norm are the average length of and statistical error associated 
with the MSTs of the random set, respectively. If /massive is significantly shorter than (I norm ) (i.e., 
by more than <r norm ), the massive stars are more concentrated and, therefore, the cluster is mass 
segregated. 

However, it is not mathematically convenient to determine the average A from a set of simula- 
tions, b ecause /massive find (/norm) from different samples may be characterized by different weights 
when averaging the value of f norm) (i.e., A). That is, a sample with very large A may dominate the 

'massive 

results. One reasonable approach is to average log(A) values instead of A. 

Therefore, we use a new variable (A') to quantify the degree of mass segregation, 

./ , /\'norm/\ , , / ^norm \ /.n 

A =log ±log . (4) 



/massive ' ^ /mi 



We will show the benefits of using the new variable in §3. 



3. Analysis and results 



A cluster loses its member stars through energy exchange and tidal effects (jConverse &l Stahler 

201CJ). In general, we are concerned about a gravitationally consistent sample, that is, we cannot 



compare a simulated cluster with several stars at great clustercentric distances to a real cluster. 
Therefore, physically speaking, it is necessary to remove the stars that do not belong to the cluster. 
We will show that this effect is not negligible. 
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The evolution of the velocity distribution is also very important for understanding the dynam- 
ical evolution of the cluster. Therefore, we will discuss the influence of the choice of the initial 
velocity distribution. In addition, we will also explore the behavior of massive stars in velocity 
space. 



3.1. Effect of escaping stars 

Removing the escaping stars from the cluster is physically necessary. On the other hand, there 
could also be some p otential escapers that may significantly affect the cluster's subsequent evolution 



(jKiipper et alj|20ld ). However, the amount of mass segregation is a property of a cluster at a fixed 
time. Therefore, it is necessary to calculate the amount of mass segregation after removing the 
escapers but all member stars (no matter whether they are potential escapers) should be included 
at any given time step. From this point of view, we only need to truncate the cluster at a specific 
radius. In the following, we will discuss the influence of our treatment of escaping stars. 

First, we compare the temporal dependence of the mass-segregation ratio without removing the 
escaping stars for both A and A'. We use average A and A' values for a given set of samples. Note 
that the average of log(A) is not equal to the logarithm of the average of A. The left-hand panel 
of Fig. Q] shows an increasing degree of mass segregation of the cluster with time from a very early 
stage. We also note that the dispersion of A increases with time, which suggests that some As may 
be extremely large, while others can be less than unity (i.e., the massive stars are less concentrated 
than the randomly selected sample stars). Therefore, if A has a large dispersion, this may lead 
to an incorrect result, as we confirmed using A'. There is a turnover point at approximately 3 
Myr, shown in the right-hand panel of Fig. [TJ That is, the cluster becomes more and more mass 
segregated until an age of 3 Myr and then the degree decreases. This turnover point may be related 
to our treatment of the escapers from the cluster. After truncating the cluster at a radius of 2.5 
pc (5r v ), we analyze the massive stars and the total number of member stars in the cluster. Fig. 
[2] shows the average mass of the 10 and 20 most massive stars and the total number left in the 
cluster as function of time. The average mass of the massive stars begins to decrease almost at the 
same time as the turnover point develops, while the cluster starts to lose members a little (< 1 
Myr) earlier. When some of the massive stars leave the cluster, the MST length that includes these 
escapers becomes very large, leading to a relatively small value of A'. This is the reason for the 
turning point in the A'-t curve. Therefore, we will radially truncate our clusters in the remainder 
of this paper. 



3.2. Choice of velocity distribution 



Because clumpy clusters evolve rapidly, small changes may significantly affect the results on 
short timescales. We find that the A' profile (as a function of time) changes if we vary the velocity 
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Fig. 1. — Evolution of the mass-segregation ratio for the 10 and 20 most massive stars (Nmst) using 
A (left) and A' (right). The line A = 1 (or A' = 0) represents the no-mass-segregation situation. 




Fig. 2. — Evolution of the average mass of the 10 and 20 most massive stars (left) and the total 
number of stars in the cluster after truncating the cluster at a radius of 5r v (2.5 pc; right). 
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distribution. In Fig. [3j we compare the results of the two different initial velocity distributions from 
Eq. (2). One is uniform throughout, while the other is radius dependent. For the radius-dependent 
velocity distribution, the degree of mass segregation reaches a peak at 0.5-1 Myr. Shortly thereafter, 
it drops to a relatively low level (similar to that of a uniform distribution). We will now attempt 
to offer an explanation for this behavior. 

We define the mean separation between stars using the MST length, 

/« MST 

= 7> 5 

n — 1 

where n is the number of stars in the sample of interest. The mean separation of a sample of 
randomly selected stars in the cluster can be a trace of the mean separation within the cluster. Fig. 
U] shows that the mean separation in the cluster and of the 10 most massive stars varies with time. 
It seems that for the radius-dependent velocity distribution, the curve of the mean separation of 
massive stars forms a valley at approximately 0.5 Myr, while the curve of the mean separation of 
the full cluster increases from the onset. This valley corresponds to core collapse of (only) the most 
massive stars in the cluster. A series of radii of different mass shells (increasing fractions of 'most 
massive' stars, from bottom to top) is presented in the left-hand panels of Fig. The right-hand 
panels of Fig. [5] show increasing fractions of 'inner cluster' stars from bottom to top. The evolution 
of the half-mass radius indicates that the cluster is undergoing a core-collapse process (driven by 
the most massive stars) until approximately 0.5 Myr, when it begins to re-expand. We compare 
this timescale to the crossing time, t croS sj which is simply the cluster radius divided by its velocity 
dispersion. For our N = 1000, D = 2.0, and q = 0.30 initial conditions, t cr0 ss ~ 0.2-0.6 Myr for the 
cluster as a whole. That is, core collapse of the few most massive stars happens during the first 
few crossing times. However, the crossing time for the core region depends on the adopted form 
for the radial velocity dispersion: it is 0.1-0.3 Myr for the uniform distribution and 0.08-0.2 Myr 
for the radius-dependent velocity dispersion. 

For the radius-dependent velocity distribution, expansion of the outer parts of the cluster after 
formation can be traced from the evolution of the radius containing 90% (by number) of cluster 
stars. Therefore, the large A' seen at early times in Fig. [3] is caused by expansion of the cluster 
rather than by a concentration of the massive stars. This expansion is caused by the initially large 
velocity in the central parts. Therefore, the radius-dependent velocity dispersion would increase 
the velocity in the central part without changing the initial virial ratio. On the other hand, the 
evolution of the MST length of the cluster as a whole, which assigns all members of a sample equal 
weights, can also be a tool for measuring the size of the cluster. 



3.3. Mass segregation in velocity space 



The massive stars are becoming centrally concentrated very rapidly, thus leading to rapid 
dynamical mass segregation. On the other hand, stars of different masses may also exhibit different 
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Evolution of mass-segregation of 1 most massive stars Evolution of mass-segregation of 20 most massive stars 




time (Myr) time (Myr) 



Fig. 3. — Evolution of the (left) 10 and (right) 20 most massive stars in the cluster for uniform and 
radius-dependent velocity dispersions. 



Evolution of mean separation Evolution of mean separation 

for radius independent velocity distribution for radius dependent velocity distribution 




time (Myr) time (Myr) 



Fig. 4. — Evolution of the mean separation for (left) uniform and (right) radius-dependent velocity 
dispersions. The solid line is the mean separation in the entire cluster, while the dashed line 
represents the mean separation of the 10 most massive stars in the cluster. 
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Fig. 5. — Evolution of (left) mass and (right) number shells for clusters truncated at 5r v . The 
top and bottom panels are characterized by uniform and radius-dependent velocity dispersions, 
respectively. The numbers in the panels refer to the fractions (by number) of cluster stars included 
in the individual curves: the left-hand panels show increasing fractions of 'most massive' stars from 
bottom to top, while the right-hand panels show the equivalent behavior for increasing fractions of 
'inner cluster' stars. 
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velocity dispersions. We adopt our MST method here to measure mass segregation in velocity- 
dispersion space. The only difference is that 'distance' now refers to the separation in velocity 
dispersion of two stars. Velocity dispersion is defined as 

o"i = Vi - v, (6) 

where <Xj is the velocity dispersion of star i, Vj the velocity of star i, and v the velocity of the whole 
cluster. Therefore, the separation in velocity dispersion of two stars can be expressed as 

<Tij = <Ti - (Tj = Vi - Vj, (7) 

which is same as separation in terms of velocities. From Fig. [6] we find that clusters can exhibit 
'inverse mass segregation' in velocity space by calculating A' using the same treatment as before. 
The reason for this is that the massive stars are preferentially concentrated in the core and have a 
higher velocity dispersion than a sample of randomly selected stars. 



4. Conclusions 

We have presented simulations and analysis of the early evolution of initially cool and clumpy 
star clusters by performing a large number of iV-body simulations. An MST method is used to 
measure the degree of mass segregation. In addition, it can also be a ruler to trace cluster size. 
Removal of the ejected stars is important physically, since it has an effect on measuring the cluster's 
mass segregation. 

We also conclude that the velocity distribution has an impact on mass segregation during core 
collapse. For a radius-dependent velocity distribution, which leads to a high velocity in the inner 
regions, the picture of (the dynamical evolution of) the cluster during the core-collapse period is 
very different from that resulting from the assumption of a constant velocity dispersion. Although 
the high-mass stars still undergo core collapse, the entire cluster expands from its formation epoch. 
As a result, a high degree of mass segregation results at the time of high-mass core collapse. We 
also adopt the MST method to discuss mass segregation in velocity space, finding 'inverse mass 
segregation,' which may be caused by the high velocity dispersion of the massive stars. 

Because our aim in this paper was to explore the effects of a radially dependent velocity 
dispersion on early mass segregation, we have not included the effects of primordial binary stars 
(beyond their dynamical formation in our simulations). This is clearly an important omission, which 
we intend to address in our future w ork. A high initial b inary fraction would introduce, on average, 



a higher velocity of escaping star s (IWeidner et al.ll2010t ). which is commonly seen in observations 



of young star clusters. However, IWeidner et al.l (|2010l ) also found that the degree of mass loss is 



nearly independent of whether or not the cluster possesses primordial binaries. In addition, for 
a N = 1000 cluster, varying the binary fraction within reasonable bounds has little effect on the 
mean mass of escapees. As a result, inclusion of initial binaries will likely not significantly affect 
the signature of mass segregation in our simulated clusters. 
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Fig. 6. — Evolution of mass segregation in velocity space. The solid and dashed lines represent 
mass segregation for a uniform, random velocity dispersion and for the radius-dependent dispersion, 
respectively. 
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